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1. INTRODUCTION 

in 

In recent decades, the experimental achievements in the investigation 
t^> of ultra-cold alkali gases^^ revived the interest to theoretical models of 

bosonic systems^. The quantum field-theoretical techniques occupy a pre- 
dominant position within the methods for the studies of weakly-interacting 
bosons. The creation-annihilation operator formalism was directly applied 
in Refs. I4I5I6I The use of nonlinear Schrodinger equation (Gross-Pitaevskii 
equation) is much more popular for the studies in this domain 7 . Some 
modifications^ and extensions^ of this methods are also available. Numerical 
techniques were utilized by KrautbP^, Pearson et alP^ and may others. The 
low-dimensional Bose-systems became of special interest recently^ 11 * 12 * 13 * 14 ^. 

In this work, a method similar to that of Bogoliubov^ is developed for a 
weakly-interacting D-dimensional Bose-gas in the harmonic trapping poten- 
tial. The paper is organized as follows. In the next Section, the Hamiltonian 
of the system under study is written in the approximate second quantiza- 
tion approach. As a result of its diagonalization a set of matrix equations is 
obtained. The results of calculations of the condensate fraction and energy 
are given in Section [371 
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2. HAMILTONIAN DIAGONALIZATION 

We consider a D-dimensional system of N weakly-interacting bosons of 
mass m confined in the harmonic trap 

V(xi, ...,x D ) = — {ulx\ + . . . + lo 2 d x 2 d ) . (1) 

The potential of interatomic interaction is given by U(x\, . . . ,xjj) = gS{x), 
where the vector x = (x±, . . . , xd), g is the interaction strength. 
The Hamiltonian reads 

+ Y, U(x j -x l )=H + U. (2) 

l<j<l<N 

Here, Pj is the momentum operator of the jth particle, Xj is its coordinate. 

It is possible to develop the second quantization approach using the 
eigenfunctions |n) = |m, . . . , no) of the operator H , i.e., an ordinary D- 
dimensional harmonic oscillator. For simplicity, we further demand that no 
ratio of the frequencies ui, . . . , wc is a rational number to avoid an accidental 
degeneration of energy levels. 

Let ali,a n be the creation and annihilation operators for the state |n). 
The corresponding energy levels are e n = h(u>ini + . . . +uz>nz)) . In this 
representation the Hamiltonian ([2]) is 

H = Y £ na[ l a ri + - Y (mnlC/lm'n^a^a^am/an/. (3) 

n m,m',n,n' 

The operators satisfy standard bosonic commutation relations: 

= S nn > (4) 

Now we apply an approximate second quantization procedure following 
Bogoliubov^. Let iVo be the number of particles at the lowest energy level 
£q. As the interaction is weak, the behaviour of bosons does not differ 
much from that of an ideal system. That is, one can expect the Bose- 
Einstein condensation to occur at low temperatures. The number iVo is thus 
a macroscopic number. As it is an eigenvalue of the operator a^ao, one 
can treat dj and clq as c-numbers: ajao = No, do a J = iVo + 1 ~ Nq, dj ~ 
\/Nq, do — y/No. To obtain more general results we do not put here Nq = N. 
Note, that for D < 3 the condensation appears only in traps^l. 

Further, following Bogoliubov, we neglect the terms having more than 
two operators with non-zero index. As the eigenfunctions |n) are real the 
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matrix elements in the second item of © at different operator combinations 
are equal: (mO|C7|nO) = (mn|J7|00) = (00|J7| mn) — g c mn , where 

= frm L \D' 2 f\( D^+n^ 1 ( m j+ n j + l \ 

Cmn -\2^h) 1 = 1 ( lj V -7777T777T V. 2 J [b) 

if rrij + rij is even for all j and c mn = otherwise. The notation cu = 
..uj d )d. 

The matrix elements with three zeros equal 

i=l v j ■ 

if nj is even for all j. 

The Hamiltonian ([3]) becomes 

H = const + ^ e n «na n + ^ ^ n + ^ n ) 

n n 

+ 9 ~Y ^2 Cmn ( 4 i«n + a^a^ + a m a n ) , (7) 



m,n 



where 'const' denotes the items of a non-operator nature. Hereafter, it will 
be dropped. For brevity, the conditions n / 0, m ^ is not written explic- 
itly. Note the appearance of linear terms ~ , ~ a which were absent in 
Bogoliubov's approach^ due to the momentum conservation law. 

Let the indices m, n run over all the states denoted by vector indices 
m, n. In order to obtain an energy spectrum of the Hamiltonian from Eq. (|7|) 
we will try to diagonalize it. For this purpose, it is possible to write (|7|) in 
the following matrix form: 

H = a f £a + 2X^/I% (a f d + a T d) + A (Wca + a f Ca tT + a T Ca) (8) 

In the above equation, a and d are vectors of infinite dimension: 

a T = (oi, a 2 ■■■), a 1 " = (a\, a\, . . . ) , d T = (di, d 2 , ■ ■ ■) , 



The diagonal matrix £ 



/ ei . . 

£ 2 . . . | anc j the matrix elements of 

v '■■ 

C are the coefficients c mn . We have also written g Nq /2 = A for brevity. 
In order to obtain the Hamiltonian in a diagonal form, 

H = ^e n o4a n = atg a, (9) 



A. Rovenchak 



where <B is the diagonal matrix with elements e n , one can apply a general- 
ization of the well-known Bogoliubov's u—v transformation demanding 

a = + Ya fT + z, a f = a j X + a T Y + z T . (10) 

Here, X and Y are square matrices of infinite dimension, we will require 
them to be Hermitian (symmetric and real), and z is a vector with real 
components. From the commutation relation (|4|) one can show that matrices 
X and Y obey the following condition 

X 2 -Y 2 = I, (11) 

where I is a unit matrix. 

It is easy to show that the linear terms (the expressions in the first 
parenthesis) of ([8]) produce only a shift of the energy levels but do not affect 
the distance between them. That is why hereafter we consider only the 
quadratic in a and a terms of the Hamiltonian. 

To eliminate the linear in a terms from H, the vector z must be 

z = -2A v / iVo(£ + 6AC')~ 1 d. (12) 

Demanding the terms atn&n and a m a n to vanish in the Hamiltonian 
one gets the following matrix equations: 

XEY + AXXCY + XXCX + XYCY = 0, (13) 
YS X + AXYCX + XXCX + XYCY = 0. (14) 

The matrix <£ is 

€ = X8X + Y8Y + AX(XCX + YCY) + 2X(XCY + YCX) (15) 

and its eigenvalues define the energy spectrum. 

If the interatomic interaction is turned off (g = 0), the solutions are 
X = I, Y = 0. We will expand the matrices X and Y into series over A: 

X = I + 2AV + Y = Xv + X 2 v 1 + ... . (16) 

The matrices Xi v i v i can De found from Eqs. (|14|) : 

v = 2 X ; x = ~£~ 1 C/2; v x = ^e- x CE- x C. (17) 

In the approximation up to A 2 we obtain the matrix (£: 

A 2 

£ = £ + 4AC + y \ {S- X cf S-ZCS- 1 C-28- 1 C 2 } (18) 
and the energy levels are given by e n = e n + 4Ac nn + 0(A 2 ). 
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Equations (|lip - (|14|) are difficult to solve, being non-linear with respect 
to infinite matrices. This set may be treated as coupled algebraic (matrix) 
Riccati equations: 

XAY + XBX + YBY = 0, (19) 

YAX + XBX + YBY = 0, (20) 

X 2 - Y 2 - I = 0, (21) 

where A = £ + 4AC, B = AC. The numerical (non-perturbative) solution 
of this set needs special approaches 16 . This problem will be considered in a 
separate paper. 

3. CALCULATION RESULTS 

If Bose-Einstein condensation occurs, the chemical potential /i of the 
system approaches zero and the total number of particles is given by 

»-» < 22 » 

where Nq is the occupation of the lowest energy level. Therefore, given N 
and the energy spectrum e n = e n (A) = e n (gNo) from the eigenvalues of ([15]) . 
one can calculate iVo as a function of temperature T. We thus arrive at a 
self-consistent problem. Its solution allows one to calculate the energy 

^ exp(e n /T) - 1 

and other thermodynamic functions. 

In Fig. Q] the dependencies of the condensate fraction Nq/N and energy 
{E — Eq)/N on temperature are given for the ID gas. For simplicity, the 
following values of parameters are used in the calculations: h = u = 1, m = 
2tt 2 , N = 1000, g = 0.0002. 

Since c nn > 0, the interaction shifts the energy levels up, which effec- 
tively corresponds to lighter particles. Thus, the Bose-condensation temper- 
ature increases. Note, however, that the presented approach is not valid for 
higher temperatures as more items must be included in the Hamiltonian ([7]). 

To summarize, the Bogoliubov's idea of an approximate second quanti- 
zation approach is generalized for a weakly-interacting Bose-gas confined in 
a harmonic trap. As a result, the equations permitting to diagonalize the 
Hamiltonian and to calculate the elementary excitation spectrum are ob- 
tained. A perturbative solution is applied and the condensate fraction and 
energy dependencies on temperature are calculated for a model system to 
demonstrate the validity of the approach. 
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Fig. 1. (left) — Condensate fraction No/N as a function of temperature T; 
(right) — energy as a function of temperature. Dashed line is the interacting 
system, solid line is the ideal Bose-gas. 
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